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basis of functions is used to form a finite-dimensional process and to 

define the Wiener process as the almost sure path- wise limit of when A'^ 
tends to infinity. We generalize such a construction to the class of centered 
Gaussian Markov processes X which can be written Xt = g{t) ■ J^^ f{t) dWt 
with / and g being continuous functions. We build the finite-dimensional 
process X^ so that it gives an exact representation of the conditional 
expectation of X with respect to the filtration generated by {^j./2'' } ^ 5: 
k<2^ . Moreover, we prove that the process X^ converges in distribution 
toward X. 
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1. Introduction 

Given some probability space, it is often challenging to establish results about 
continuous adapted stochastic processes. As a matter of fact, the mere exis- 
tence of such processes can prove lengthy and technical: the direct approach 
to build continuous Markov processes consists in evaluating the desired finite- 
dimensional distributions of the process, and then constructing the measure 
associated with the process on an appropriate measurable space, so that this 



measure consistently yields the expected finite-dimensional distributions (|l3l ). 
In that respect, it is advantageous to have a discrete construction of a contin- 
uous stochastic process. For more general purpose, a discrete representation of 
a continuous process proves very useful as well. Assuming some mode of prob- 
ability convergence, at stake is to write a process X as a convergent series of 
random functions /„ • ^„ 

oo N 

Xt = y^ Ut) ■ in = lim V Ut) ■ in , 
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where /„ is a deterministic function and ^„ is a given random variable. 
The Levy construction of Brownian motion -later referred as Wiener process- 
provides us with a first example of discrete representation for a continuous 
stochastic process. Noticing the simple form of the probability density of a Brow- 
nian bridge, it is based on completing sample paths b y in terpolation according 
to the conditional probabilities of the Wiener process (jlOl ). More especially, the 
coefficients ^„ are Gaussian independent and the elements /„, called Schauder 
elements, are obtained by time-dependent integration of the Haar elements. This 
latter point is of relevance since, for being a Hilbert system, the introduction of 
the Haar basis greatly simplify the demonstration of the existence of the Wiener 
process (0). 

From another perspective, fundamental among discrete representations is the 
Karhunuen-Loeve decomposition. Instead of yielding a convenient construction 
scheme, it represents a stochastic process by expanding it on a basis of or- 
thogonal functions (0; [ll|). The definition of the basis elements /„ depends 
only on the second-order statistics of the considered process and the coefiicients 
are pairwise uncorrelated random variables. Incidentally, such a decompo- 
sition is especially suited to study Gaussian processes because the coefiicients 
of the representation then become Gaussian independent. For these reasons, 
the Karhunen-Loeve decomposition is of primary importance in exploratory 
data analysis, leading to methods referred as "principal component analysis", 
"Hotelling transform" (@) or "proper orthogonal decomposition" |T3 ) according 
to the field of application. In particular, it was directly applied to the study of 
stationary Gaussian Markov processes in the theory of random noise in radio 
receivers (0). 

In view of this, we propose a construction of Gaussian Markov processes us- 
ing a Haar-like basis of functions. The class of processes we consider is general 
enough to encompass commonly studied centered Gaussian Markov processes 
that satisfy minimal properties of continuity. We stress the connection with the 
Haar basis because our basis of decomposition is the exact analog of the Haar- 
derived Schauder functions used in the Levy construction of the Wiener process. 
As opposed to the Karhunene-Loeve decomposition, our basis is not made of 
orthogonal functions but the elements are such that the random coefficients ^„ 
are always independent and Gaussian with law A/'(0, 1), i.e. with zero mean and 
unitary variance. 

The almost sure path-wise convergence of our decomposition toward a well- 
defined continuous process is quite straightforward. Most of the work lies in 
proving that the candidate process provides us with an exact representation of 
a Gaussian Markov process and in demonstrating that our decomposition con- 
verges in distribution toward this representation. Validating the decomposition 
essentially consists in proving the weak convergence of all the finite-dimensional 
measures induced by our construction on the Wiener space: it requires the in- 
troduction of an auxiliary orthonormal system of functions in view of using the 
Parseval relation. To furthermore establish the convergence in distribution of 
the representation, we only need demonstrating the tightness of this family of 
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induced measures. 

The discrete construction we present displays both analytical and numerical 
interests for further applications. 

Analytically- wise, even if it does not exhibit the same orthogonal properties 
as the Karhunene-Loeve decomposition, our representation can prove as ad- 
vantageous to establish analytical results about Gaussian Markov processes. It 
is especially noticeable when computing quantities such as the characteristic 
functional of random processes ([g; y) as shown in annex. This is just an ex- 
ample of how, equipped with a discrete representation, one can expect to make 
demonstration of properties about continuous Gaussian Markov processes more 
tractable. If the measures induced by our decomposition on the classical Wiener 
space converge weakly toward the measure of a Gaussian Markov process, we 
put forward that the convergence of our decomposition is almost sure path- wise 
toward the representation of a Gaussian Markov process. This result contrasts 
with the convergence in mean of the Karhunene-Loeve decomposition. 
From another point of view, three Haar-like properties make our decomposition 
particularly suitable for certain numerical computations: all basis elements have 
compact support on an open interval with dyadic rational endpoints; these in- 
tervals are nested and become smaller for larger indices of the basis element, 
and for any dyadic rational, only a finite number of basis elements is nonzero 
at that number. Thus the expansion in our basis, when evaluated at a dyadic 
rational, terminates in a finite number of steps. These properties suggest an 
exact schema to simulate sample paths of a Gaussian Markov process X in an 
iterative "top-down" fashion. Assuming conditional knowledge of a sample path 
on the dyadic points of Dn = {fc2~^|0 < k < 2^}, one can decide to further 
the simulation of this sample path at any time t in Dat+i by drawing a point 
according to the conditional law of Xt knowing {ATtj^g^^, which is simply ex- 
pressed in the framework of our construction. It can be used to great advantage 
in numerical computations such as dychotomic search algorithms for first pas- 
sage times: considering a continuous boundary, we shall present elsewhere a fast 
Monte-Carlo algorithm that simulates sample-paths with increasing accuracy 
only in time regions where a first passage is likely to occur. 



2. Main Result 

Beforehand, we emphasize that the analytical and numerical advantages granted 
by the use of our decomposition come at the price of generality, being only 
suited for Gaussian Markov processes with minimal properties of continuity. We 
also remark that if the Karhunen-Loeve decomposition is widely used in data 
analysis, our decomposition mainly provides us with a discrete construction 
scheme for Gaussian Markov processes. 

Proposition. Let X = {A't, J^t;0 <t <l} be a real adapted process on some 
probability space {Q,J-,P) which takes value in the set of real numbers and let 
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^n,k with n>0 and 0<fc<2" be Gaussian random variables of law Af{0, 1) . 
If there exist some non-zero continuous functions f and g such that 

Xt = g{t)-[ f{t)dWt, with 0<t<l, 
Jo 

tlicii there exists a basis of continuous functions ^n,k for n>0 and 0< A;<2" 
such that the random variable 

N 

n=0 0<A;<2'»-i 

follows the same law as the conditional expectation of Xt with respect to the 
filtration generated by for 0<k<2^ . The functions 'i'n.k thus defined 

have support in Sn,k = (fc + l)2^"+"'"] and admit simple analytical ex- 

pressions in terms of functions g and f. 

Moreover, the path-wise limit liniAr^^oo -''^^ defines almost surely a continuous 
process which is an exact representation of X and we have 

meaning that the finite- dimensional process X^ converges in distribution toward 

X when N tends to infinity. 

Remark. The function f can possibly be zero on a negligible set in [0, 1] in the 
previous proposition. 

To prove this proposition, the paper is organized as follows. We first review 
some background about Gaussian Markov processes X and their Doob represen- 
tations as Xt = g{t)- Jq f{t) dWt- Then we develop the rationale of our construc- 
tion by focusing on the conditional expectations of the process X with respect 
to the filtration generated by for < fc < 2^. In the fifth section, we 

propose a basis of expansion to form the finite-dimensional candidate processes 
X^ and justify the limit process X as the almost sure path-wise convergent 
process YanM^oo X^ . In the sixth section, we introduce the auxiliary Hilbert 
system and prove an important intermediate result. In the last section, we show 
that the finite-dimensional processes X^ converge in distribution toward X and 
that X is an exact representation of X . 

3. Background on Gaussian Markov Processes 
3.1. Basic Definitions 

We first define the class of Gaussian Markov processes. Let us consider on some 

probability space (fJ, P) a real adapted process X = {Xi,Ti \ <f <oo} which 
takes value in the set of real numbers. We stress that the index t of the random 
variable Xt runs in the continuous set For a given realization a; in fi, the 
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collection of outcomes t Xt{uj) is a. sample path of the process X. We only 
consider processes X for which the sample paths 1 1-^ Xt{uj) are continuous. With 
these definitions, we are in a position to state the two properties characterizing 
a Gaussian Markov process. 

1. We say that X is a Gaussian process if, for any integers k and positive 
reals ti < t2 < ■ ■ ■ < tk, the random vector {Xt-^ , Xt^ , • ■ • , Xt^ ) has a joint 
normal distribution. 

2. We say that X is a Markov process if, for any s,t>0 and F e S (M), with 
B (M) the set of real Borelians, 

P iXt+, eT\T,)^P {Xt+s e r I X,) , 

which states that the conditional probability distribution of future states 
Xt+s, given the present state and all past states J-'g, depends only upon 
the present state Xs- 

A Gaussian Markov process is a stochastic process that satisfies both Gaussian 
and Markov properties. 

The Wiener process and the Ornstein-Uhlenbeck process are two well-known 
examples of Gaussian Markov processes. The Wiener process is defined as the 
only continuous process Wt for which Wq = and the increments Wt — Wg 
are independent of Ts and normally distributed with law A/'(0,t — s). These 
requirements naturally place the Wiener process in the class of Gaussian Markov 
process. The Ornstein-Ulhenbeck process can be defined as a solution of the 
stochastic differential equation of the form 

dXt = aXt dt + dWt with a e M , 

We designate C/" the Ornstein-Ulhenbeck process of parameter a starting at 
for i = and we give its integral expression 

C/f = / e"("-*) dWu . (1) 



The process t/" naturally appears as a Gaussian Markov process as well: it is 
Gaussian for integrating independent Gaussian contributions and Markov for 
being solution of a first-order stochastic differential equation. It is known that 
both processes can be described as discrete processes with an appropriate basis 



of random functions (|14| ) . 



3.2. The Doob Representation 

The discrete construction of the Wiener process and the Ornstein-Uhlenbeck 
process is likely to be generalized to a wider class of Gaussian Markov processes 
because any element of this class can be represented in terms of the Wiener 
process. By Doob's theorem for any Gaussian Markov process X, there 
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exist a real non-zero function g and a real function / in Lf^^{R^) such that we 
have the integral representation of X 

Xt^g{t)- I f{t)dWt. (2) 

where W is the standard Wiener process. If we introduce the non-decreasing 
function h defined as ^ 

h{t) = f f{u)du, 
Jo 

then, for any t, s > 0, the covariance of X can be expressed in terms of functions 
h and g 

E {Xt ■ X,) = g[t)g{s) ■ h{ min(t, s)) . (3) 

The Doob's representation ^ indicates that X is obtained from by a change 
of variable in time t ^—^ h(t) and, at any time t, by a change of variable in 
space by a time-dependent factor x ^ g{t) ■ x. The couple of functions (/, g) 
that intervenes in the Doob's representation of X is not determined univocally. 
Yet, one can defined a canonical class of functions (/, g) which are uniquely 
defined almost surely in Lf^^iW^) if we omit their signs (fsl). Incidentally, we 
can compare the integral formulation ([2]) with expression ^ of the Ornstein- 
Uhlenbeck process: we remark that the representation of this Gaussian Markov 
process is provided by setting g{t) = e"* and f{t) = e~"*, which happens to be 
its canonical representation. 



3.3. Analytical Results 

The discrete construction of Gaussian Markov processes will rely on two ana- 
lytical results that we detail in the following. 

First, the Doob's representation allows us to give an analytical expression for 
the transition 'keinelp{Xt=x \Xt„ = xo) of a general Gaussian Markov process. As 
the Doob's representation is a simple change of variables, it is easy to transform 
the expression of the Wiener transition kernel p {Wt=x \ Wtg = xo) to establish 



p{Xt=x \Xta = XQ) = , • exp 



2 - 



We have to mention that this expression is only valid if h{t) ^ h{to), otherwise 
X is deterministic and p{Xt=x \Xtg = xo) = 5xo{x). 

Second, we can use this result to evaluate q{Xt^ = y \Xt^= x, Xt= z) with < 
ty<tz, the probability density of Xf knowing its values x and z at two framing 
times and t^ - Because X is a Markov process, a sample path 1 1-^ Xt{oj) which 
originates from x and joins z through y is just the junction of two independent 
paths: a path originating in x going to y and a path originating from y going 
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to z. Therefore, after normalization by the absolute probability for a path to go 
from X to y, we have the probability density 

fv w V \ P{Xty=y\Xt=x)-p{Xt=z\Xt=y) 
qiX,^^y\X,^^x,X,^^z) = piX,^.^\X,^.x) ' 

Thanks to the previous expression, we can compute the distribution of Xt^ 
knowing Xt^ and Xt^ , which is expected to be a normal law because we only 
consider Gaussian processes. For a general Gaussian Markov process X, we 
refer to that probability law as A/'(xM(^i;): with mean value x^^{ty) 

and variance x<^{tyY- We show in annex that these parameters satisfy 

„/, N ^ gftjy) - Hty) I ajty) K^y) - K^x) 

git.) ' hit,) ~ h{t,) ' ^ g{t,) ' h{t,) ~ h{t,) ' ' ' 

I, ,2 2(, . {hity)-hit,)){h{t,)-hity)) 

= ^ ^'y^ Mt.)-Mt.) ■ ^ ^ 

Once more, we have to mention that these expressions are only valid if hity) ^ 
hitx). Wether considering a Wiener process or an Ornstein-Uhlenbeck process, 
the evaluation of (|4]) and (O with the corresponding expression of g and h leads 
to the already known results 



4. The Rationale of the Construction 



4-1. Form of the Discrete Representation 



Form now on, we will suppose that the zeros of the function / pertain to a 
negligible ensemble in [0,1], causing the function h to be strictly increasing. 
We will further restrain ourselves to Gaussian Markov processes for which the 
functions / and g belong to the set of continuous functions on [0, 1] denoted 
C(0, 1). We remark that in such a case, the functions 1 1— s- xfJ'it) and 1 1-^ xf (t) 
are continuous on [0, 1]. 

Bearing in mind the example of the Levy construction for the Wiener process, 
we want to define a basis of continuous functions 'i'n.k in ^(0, 1) with < A; < 2" 
to form the discrete process 

N 

^t'^^E E '^nAt)-^n,k, 
n=0 0</c<2i-i 

where £^n,k are independent Gaussian random variables of standard normal law 
A/'(0, 1). We want to chose so that t i-^ X^iuj) converges almost surely 

toward 1 1—>- Xiuj) when N tends to infinity. Given the continuous nature of the 
processes X^, we require that the convergence is uniform and normal on [0, 1] 
to ensure the definition of a continuous limit process X = limAr^oo X^ on [0, 1]. 
Moreover, we want to define '^n,k on supports Sn.k of the form 



5„,fc=[fc.2-"+i,(fc + l)2-"+i] 
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As a consequence, the basis of functions "^n.k will have the following proper- 
ties : all basis elements have compact support on an open interval with dyadic 
endpoints; these intervals are nested and becomes smaller for larger indices n 
of the basis element, and for any dyadic rational, only a finite number of basis 
elements is nonzero at that number. 

4-2. Conditional Averages of the Process 

Now remains to propose an analytical expression for '^n,k- If we denote 
the set of reals {fc2~^ | < k < 2^}, the key point is to consider = 
E{Xt\{Xs}, s £ Dn) the conditional expectation of the random variable Xt 
given Xs with s pertaining to the set of dyadic points D^. The collection of 
random variables defined on fl specify a continuous random process Z^ on 
fi. We notice that, if t^ = fc2"^ and t,_ = {k+l)2-^ with 0<fc<2-^ are the two 
successive points of Dj^ framing t, the random variable Zj^ is only conditioned 
by Xt^ and Xt^ : 



Using expression (j4|) , we can express the sample paths t Z^ {us) as a function 
of t on [tx, tz]'. for a given uj in the sample space fl, we write 



where the conditional dependency upon parameters Xt^{uj) = x and Xt^iw) — z 
is implicit in xl^-^'^ ■ The random process Z^ appears then as a parametric func- 
tion of {Xs\ g^jj^'. for any uj in the sample space $7, the sample path t ^ Xtito) 
determines a set of value {xs}^,^^,^ = {Xgiuj)} ^^j^^ and by extension a sample 
path t ^ Z^ {uj) for the process Z^ . 

Now two points are worth noticing: first, t ^ Xt{ijj) and t ^ Zf{uj) are contin- 
uous sample paths that coincide on the set D^; second, we have {0, 1} = Dq C 
-Di C • • • C Dat a growing sequence of sets with limit ensemble V the set of 
dyadic points in [0, 1], which is dense in [0, 1]. Then, provided the path-wise con- 
vergence is almost surely normal and uniform on [0, 1], the limit process of Z^ 
when N tends to infinity should be continuous and the processes limAr^oo Z'^ 
and X should be indistinguishable on f2. 

4-3. Identification of Conditional Averages and Partial Sums 

Identifying the process Z^ with the partial sums X'^ provides us with a ratio- 
nale to build the functions ^n,k- 

We first need to consider the random variable Z^^'^ on the support Sn+im = 
[tx,tz] of the function "i'pf^i^k- The Markov property of the process X entails 



■N 



E{Xt\{X,},se Dn) = E{Xt\Xt^,XtJ. 




. X def N k / 1\ 
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Hence, for a given u the estimation of the sample paths t Z^~^^ (uj) is now 
dependent upon y = Xty{u)) with ty the midpoint of tx and 

,N+l,,,-,_fxfM^,Ly{t,X,y) if tx<t<ty def ^^N,kt 



\ xl^ty,tAt,y,Z) if ty<t<t^ ^'^^ 

Wc can identify the conditional process and the partial sums . Then, for 
any lu in Q, writing the sample path 1 1— > Z^^^{lu) as a function of y, we have 

Assuming conditional knowledge on Dn, the quantity Z^ =E{Xt\{Xs}, s G D^) 
becomes deterministic and the outcome of the random variable is only 

dependent upon the values of the process on Dat+i \ D^. More precisely, on the 
support SN+i,k, the outcome of Z^^^ is determined through the function xV^''^ 
by y the outcome of Xt^ given Xt^=x and Xi^=z. 

The distribution of X^^^ given Xt^ = x and Xt^ = z follows the law M{xl^{iy), 
xO'{ty)) and we denote Ypf^k a Gaussian variable distributed according to such a 
law. With this notation, we are in a position to propose the following criterion 
to compute the function ^'jv+i,fe: the element n+im is the only positive function 
with support included in SN+i,k such that the random variable N+i,k(t) ■ £.N+i,k 
has the same law as ^Ar,fc) — x^J'^'''{t)■ Direct calculations confirms 

that the previous relation provides us with a consistent paradigm to define 
the functions \l/jv_)_i^/j. Incidentally, we have an interpretation for the statistical 
contribution of the components ^'iv+i,fc • ^N+i,k' if one has previous knowledge 
of Xt on Dat, the fimction ^N+i.k ■ (.N+i.k represents the uncertainty about 
Xt that is discarded by the knowledge of its value on -D/v+i \ -D/v. 

5. The Candidate Discrete Process 
5.1. The Basis of Functions 

We recall that we carry out the case for which the function / has a negligible 
set of zeros in [0,1], which directly follows from the previous section. Before 
specifying the candidate basis elements we introduce the following short 

notations to simplify the writing of their expressions 

«„,fc = (2/e) 2-" , m„,fc = (2A:+l)2-" , r„,fe = 2 (fc+1) 2"" . 

Then for n > and < fc < 2"~^, the explicit formulation of the basis of 
functions k reads 



^n,k{t) 



Ln,k ■ g{t){h{t) - h{ln,k)) if ln,k <t < mn,k , 

Rn,k ■ 9{t){h{rn,k) - h{t)) if mn,k<t<rn,k, (6) 
otherwise , 
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where we use the constants Ln^k and Rn.k that are defined by the relations 




h{rn,k) ~ Hrrin^k) 

(Hrn^k) - h{ln,k)) {Hmn,k) - Wn,k)) 



h{mn,k) - h{ln,k) 

{Hryuk) ~ hiln,k)){h{rn^k) ~ hi'mn,k)) 



For = 0, the basis element ^o.o needs to satisfy the relation 
*o,o(t) • eo,o = E{Xt\{X,}, seDo^ {0, 1}) , 
which completely defines the analytical expression of ^'o.o as follows 

m ■ {h{t) - fe(/o,o)) 

\/h{rQ^o) - h{lQ.o) 

As expected, we directly ascertain the continuity of the '^n,k by continuity of / 
and g. 

We should briefly discuss the form of the functions 'i'n.k- In the case for which 
g{t) = 1 and h{t) — t, we find the usual expression of 'i'n,k for the Levy construc- 
tion of a Wiener process: the elements of the basis are the triangular wedged- 
functions obtained by integration of Hn^k the standard Haar functions. In the 
general case of a Gaussian Markov process, the expression of 'i'n,k can be de- 
rived from the Wiener process basis elements by three operations: a change of 
variable in time dt i-^ h'{t)dt = p{t)dt, a time-dependent change of variable in 
space X ^ g{t) ■ x and a multiplication by the coefficients Ln,k and Rn,k- The 
effect of this multiplication by Ln,fc and Rn.k will be explain in section [7l 
Moreover, the paradigm of the construction makes no assumption about the 
form of the binary tree of nested compact supports Sn,k- Let us consider a given 
segment /q^o = [^o,Oj''o,o[ and construct by recurrence such a tree. We suppose 
that we have the following partition 

-^0,0= [J In.k— [J [^n,fe, ?'n,fe [ ■ 
0<fe<2"-i 0<fc<2"-i 

For each k such that < fc < 2"~^, we draw a point m„.fc in Sn.k- Then, we have 

In,k ~ [^n,fc,^n,fe[ [J \jnn,k,'rn,k[ = In+l,2k [J In+l,2k+l 

and by construction, we posit mn,k = ln+i,2k+i = Tn+i,2k- Iterating the process 
for increasing n, we build a tree of nested compact supports Sn,k = In,k- 
The definition ([6|) enables us to explicit elements 'i'n.k that are adapted to any 
such tree. The so-defined functions 'i'n.k will appear to be valid basis elements 
to build a discrete representation of X under the only requirement that 

lim sup \rn,k - In.kl = ■ (7) 
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5.2. The almost sure Normal and Uniform Convergence 

We want to prove the validity of the discrete representation of a Gaussian 
Markov process with Doob's representation ([2]) using the proposed basis of func- 
tions 'i>n,k- Let us consider the partial sums defined on by 

N 

^^ = Y1 J2 *n,fcW-en,fc for i e 5*0,0 - [0, 1] . 
71=0 0<fc<2"-i 

We need to study the path-wise convergence of the partial sums on fl to 
see in which sense we can consider limjv^oo as a proper stochastic process. 
Again, we only consider Gaussian Markov processes for which the functions / 
and g belong to the set of continuous functions C[0, 1]. If we designate the L°° 
norms of / and g on [0, 1] by \\f\\^ and we can show that 



sup sup < 

0<fe<2"-i 0<i<l 

< 

For (/, g)-bounded Gaussian Markov processes, this inequality provides us with 
the same upper bound to the elements ^n,k as in the case of a Wiener process 
times a constant ||5||^ • ^^"^ same Borel-Cantclli argument as for 

the Haar construction of the Wiener process (8), for almost every w in fi, the 
sample path converges almost surely normally and uniformly in t to a function 
t H-> Xt{oj) when N goes to infinity. 

It is worth noticing that, since / and g are continuous functions, so are the 
basis functions ^'„.fe- Then, for every cu in fi, the sample path t i— > X^ [uj) is 
a continuous function in C[0, 1]. As the convergence when N tends to infinity 
is normal and uniform in t, the limit functions t ^ Xt{uj) results to be in 
C[0, 1] almost surely on il. This allows us to define on 17 a limit process X = 
limjv_+oo X-^ with continuous paths. 

Showing that X is an admissible discrete representation of the Gaussian Markov 
process X only amounts to demonstrate that, for any integers k and positive 
reals ii < t2 < • • • < ife, the random vector {Xti,Xt2^ - ■ ■ ,Xt,.) has a the 
same joint distribution as {Xt-^, Xt2, ■ ■ ■ ,Xt,.). As X is defined as the path- 
wise almost sure limit of X^ when N tends to infinity, this result is implied 
by the convergence in distribution of the continuous processes X^ toward their 
limit X (@;[l!). We will therefore establish the convergence in distribution of our 
representation in section[5]and incidentally validate X as an exact representation 
of AT. In that perspective, we devote the following section to set out the meaning 
of the convergence in distribution of our candidate process X^ . 



— ■ ; — — r~ ~ ■ ■ "-^""o 

' n.k ''71, k 

2-"*'-||5lL-|l/lloo- (8) 
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6. The Convergence in Distribution 
6.1. The Finite- dimensional Measures 

In this section, we speciiy the finite-dimensional probabiHty measures in- 
duced by the processes . Beforehand, we introduce the notations 



[n,fc]=2""^ + fc with 



[0, 0] 
[0, 1] = 2^ 



to allow us to list the midpoints nin^k of the tree of supports Sn^k in the prefix 
order. By reindexing according to 

t[nM = mn^k = (2fc+l) 2"" , 

we get an ordered sequence to < ti < ■ ■ ■ < ■ Let us now consider the 
finite-dimensional space of admissible functions for 



{*n,fc}os 



C Vect I , I o<n<N 



When it is equipped with the L°° norm, is a complete, separable metric 
space under the distance d{f,g) = \\f — g\\^- We can provide the space 
with the (T-algebra B{C^) generated by the cylinder sets Cbq,- - ,b,n : which are 
defined for any collection of Borel sets Bq, Bi, - ■ ■ , B2N in B{M.) by 



C. 



{x e I x(t[„,fc]) e B[„,fc], 0<n<A^, 0<fc< 2"-i} 



The random process X^ induces a natural measure on {C^ , B{C^)), such 
that 

VBeS(C^), P^{B)^P{Lu\X^{Lu)eB). 

Since for any x in we have x{0) — 0, the induced measure is entirely 
determined on the cylinder sets of the form Cbq,- - ,b^h with — {0}- Keeping 
this in mind, we show in annex that P^ admits a probability density : for 
any cyhnder set Cbo, -- ,b « of B{C^) with Bq = {0} we have 



J Bi ^ B^^ 



{xi, - • • , X2n) dxi ■ ■ ■ dx2N 



where p^ is made explicit with the help of the transition kernel p oi X defined 
in dH) 

2"-l 

p^{xi, - ■■ ,X2n) = W p{xk+l,tk+l\xk,tk) ■ 
k=0 

We want to specify in which sense the finite-dimensional probability measures 
defined on (C^,S(C^)) converge to a Hmit measure P associated to X. 
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6.2. The Weak Convergence 

Here, we consider that the stochastic processes and X take value in the 
Wiener space, that is the space of continuous functions C[0, 1]. This allows us 
to characterize the X-induccd measure fx = P associated with X on C[0, 1]. 
Defining the X^-induced measures /x^ on C[0, 1] as well, we then state the 
convergence of toward P in terms of weak convergence of /z^ toward /x. 
The process X^ defined on some probability space (fi, !F, P) have continuous 
sample paths t X^{lu) and so does the Gaussian Markov process X. Being 
a complete, separable metric space under the distance d{f,g) = \\f — g\\^, the 
Wiener space is a natural space to define X^- and X-induced measures. We 
consider X^ and X as a random variables with values in the measurable space 
(C[0,1],B(C[0, 1])), where B(C[0,1]) is the a-field generated by the cylinder 
sets of C[0, 1]: as previously X^ and X naturally induce the probability mea- 
sures and /i defined by 

li^{B) = P{ojen\X^{oj)eB} and ii{B) = P{oj & n\X{oj) & B} 

for any B in S(C[0.1]). 

More specially, the measure /x is called the Wiener measure of the Gaussian 
Markov process X. Assuming a general cylinder set to be 

^b;;::.'*s„ = e c[o, i] | xitk) e Bk, o<k<n} 

for any 0<ti<t2<- ■ ■<tn<l and any Borel sets Bi,B2, ■ ■ ■ ,Bn in B{R), is 
the unique probability measure such that for any set 

B„) = / ■ ■ ■ / p{Xl,ti\0,0) ■ ■ ■p{Xn,tn\x„-i,tn-l) dXt^ ■ ■ ■ dXt„ . 

Through the induced measures fi^ and fi, we want to study the convergence of 
the process X^ toward X from a probabilistic point of view. In that respect, 
the most general form of convergence one might expect is the weak convergence. 
We recall that /j,^ is weakly convergent to /x if and only if for any bounded 
continuous function of C[0, 1], we have 

lim / (j)(x) dfj,^ {x) = / (j){x) diJ,{x) , 
Jc Jc 

where C is a short notation for C[0, 1]. 

6.3. The Convergence in Distribution 

We can now translate on the Wiener space the much desirable property that the 
continuous processes X^ converges in distribution toward the Gaussian Markov 
process X, which we denote 
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We recall that, by definition, X converges in distribution to X if and only if 
for any bounded continuous function (j) in C[0, 1] we have 



where Epf and E are the expectations with respect to P and P respectively. 
With the definitions made in the previous section, the convergence in distribu- 
tion of the representation X^ is rigorously equivalent to the weak convergence 
of the measures /i^ toward /i on the Wiener space (C[0, 1],B (C[0, 1])). 
We will show this point in section [8] following the usual two steps reasoning 
inspired by the Prohorov theorem (1). To show the weak convergence of the 
sequence of measure fi^ , it is enough to prove the two statements: 

1. For every integer fc > and reals < ti < t2 < ■ ■ ■ < tk < 1, the finite- 
dimensional vector {X^^, X^, • • • , X^) converges in distribution to {Xt^ , 
Xt2 , • • • , Xt^. ) when N tends to infinity. 

2. The family of induced measures fj,^ is tight: for every rj > 0, there exist a 
compact K C C[0, 1] such that fj,'^ (K) > 1 — 77 for every N. 

Before establishing these two criteria, we first need to compute the limit of the 
covariance of when N tends to infinity. This calculation is the crucial point 
to validate our representation and we will carefully detail it in the following 
section. 

7. The Covariance Calculation 

7.1. Definition of the Auxiliary Basis 

Let us remember that the Gaussian Markov process X admits a Doob's rep- 
resentation ([2]). If we posit appropriate regularity properties for / and g, it 
is straigtforward to see that, by Ito formula, such a process is solution of the 
stochastic differential equation 



It is then tempting to inject the proposed basis element in the previous 

equation and to consider the functions $n,fc defined as 



hm En{cI>{X'')) ^E{^{X)) 





The functions ^n,k are actually well-defined despite the division by /, which is 
potentially zero, since calculations show that the ^n,k are given explicitly for 
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n > by 



Ln,k ■ fit) if ln,k <t< rUn^k , 

Rn.k ■ f{t) if rrin.k < t < rn,k , 
otherwise . 

As for the element ^o,o, it gives rise to the well-defineci function 

fit) 



*o,o(i) 



VK^ofi) - Kiofi) 



(9) 



We will show that the family of functions ^n,k is a Hilbert system of a subspace 
of L^(0, 1), a property that will prove useful to compute, for any t and s in [0, 1], 
the limit of the covariance EiX^ ■ X^) when TV tends to infinity. The consider- 
ation of this result will also enable us to interpret the coefficients L„.fc and Rn,k- 



7.2. Characterization as an Hilbert System 

From now on, we denote if,g) the usual inner product of / and g in L^(0, 1). 
Let us introduce £f, the closure of the vectorial space £f defined as 

- {0 e L\Q, 1) 1 e L\0, 1), = (^1 . 

With the usual inner product, the space £f inherits the structure of a Hilbert 
space from L^(0, 1) for being a closed subspace of a Hilbert space. It is immediate 
to see that the $„,fc belong to £f when written as ([9]). We need to prove that 
the '^n,k constitute an orthonormal family for the usual inner product and that 
the vectorial space of their finite linear combinations is dense in 
Let us start with the orthonormal property of the family and consider two 
functions and with (n, fc) ^ (n',/c'). If n = n' and k ^ k', ^n,k 

and ^n',k' have disjoint supports and their inner product is necessarily zero. 
Assuming that n' > n, ^n,k and ^n',k' have intersecting supports Sn,k and 
Sn',k' if and only if Sn,k is strictly included in Sn',k'- Then, it is very useful to 
remark the nullity of the following inner product 

rnin.k f^Ti.k 
if,^n,k) = Ln,k- f^iu)du + Rn^k- f^iu)du 

•lln.k JrUn.k 

_ / ihirn.k) - hiln,k)) iKmn^k) - /i(^n,fc)) 

Y hirn,k) - hiinn,k) 

_ / (fe(rn,fc) - hiln,k)) {hjmn^k) - hjln^k)) _ Q 
y hirn,k) - himn,k) 
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which entails that ^n,k and ^n',k' are orthogonal for n ^ n' . As for the norm 
of the functions ^n,k, we directly compute for n > 

{^n,k.^nM) = L^k- fiu)du + Rl,- f{u)du 

_ hjrn^k) - Krrin.k) ^ h{mn,k) - h{ln,k) _ ^ 

h{Tn,k) - Wn,k) h(rn,k) - Wn,k) 

and it is straightforward to see that (^o.o, 5'o,o) = 1 for n = 0. Hence, we have 
proved that the collection of $„ ^ forms an orthonormal family of functions in 

We still have to show that the linear combinations of ^n,k generate a dense 
vectorial space in £/. We can easily be convinced of this point once we consider 
the family of functions / • Hn^k, where the Hn^k designate the Haar functions 
adapted to the supports Sn,k - The orthonormal system Hn,k is dense in L^(0, 1) 
as soon as condition ([7|) is satisfied. Then, as each / • Hn,k can be obtained by 
finite linear combination of we conclude that 




The interpretation of the coefficients Rn,k and L„^fe is made conspicuous in this 
context. The natural surjective morphism of L^[0,1] to £f is the application 
(j) i—f f ■ (j). The family / • Hn^k in £/ is the image of the Haar basis Hn,k, but 
it is not a Hilbert system of £/. The coefficients Rn.k and Ln,k results from the 
operations of orthonormalization of the family / • Hn k to form a Hilbert system 

7.3. Application of the Parseval Relation 

Now that these preliminary remarks have been made, we can evaluate, for any t 
and s in [0, 1] , the limit of the covariance E{X^ ■ Xf) when N tends to infinity. 
As S,n,k are independent Gaussian random variables of normal law A/'(0, 1), we 
see that the covariance oi is given by 

N 

E (Xf . ) = ^ ^r.At) ■ *n,fe(s) • (10) 

n=0 0<A;<2"-i 

To compute the limit of the right-hand side in (fTO| , we need to remark that the 
element of the basis 'i>n,k and the function of the auxiliary Hilbert system ^n,k 
are linked by the following relation 

= I X[o.t]{u)g{t)f{u)-^n.k{u)du. (11) 
-'0 
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In the previous expression, we use the indicator functions of the segment [0,t] 
defined as 

I if < u < i 
otherwise 



X[o,t]i 



u 



In order to simphfy the notations, we now introduce the functions rit elements 
of the Hilbert space £/ 

Vt{u) = X[o,t](w)5W/(") ■ 

With the help of the function T]t , we can then write the integral definition of the 
basis element '^n,k (fTT|) as the inner product in £f 

*n,fe(i) = (?yt,$n,fc) • 

Remembering that the family of functions ^n,k is a Hilbert system of £f, we 
can make use of Parseval identity, which reads 

oo 

iVs,Vt)^Yl (?7s,^n,fc) • (?7t,*rKfc) . (12) 

n=0 0<fc<2"-i 

Thanks to this relation, we can conclude the evaluation of the variance of X, 
since a direct explicitation of yields 



/' 

Jo 



X[o,t]{u)g{t)f{u) ■ X[o,s]iu)g{s)f{u)du = ^ ^ *„,fe(0 ■ *n,fc(s) . 

n=0 0<fc<2"-i 

(13) 



The left term in p3|) precisely happens to be the same as the covariance of X 
given by relation ^ and we recap the statement by saying 

hm E (Xf • ) = gis)git) ■ h{ min(t, s)) ^ E {Xt ■ X,) . 



oo 



In section IHl we will use this relation to show that the random vector {X^,X^, 
■■■ ,X^) converges in distribution toward {Xt^, Xt^,--- ,Xt^) for any reals 
0<ii<t2<- • •<<fc<l. 



8. The Convergence of the Representation 
8.1. Formulation of the two Sufficient Criteria 

We are now in a position to proceed to the demonstration of the weak conver- 
gence of the induced measures /i^ toward ii. This will prove the convergence in 
distribution of the representation 
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and validate the fact that the hmit process X is an exact representation of X. We 
just need to verify that our representation satisfies the two criteria mentioned 
in section [6l 

First, we show that the family of measures iin fulfill the tightness condition, 
using a characterization of tightness on C[0, 1]. This characterization is a version 
of the Arzela-Ascoli theorem. It is formulated in term of mod{x, 5) the (5-modulus 
of continuity of a function x in C[0, 1] : 

V(5 > 0, mod(x,S) — max \x(t) ~ x(s)\ . 

\t-s\<S ' ' 

Bearing in mind that we only consider the subset of functions x in C[0, 1] for 
which a;(0) = 0, the characterization reads 

lim sup /i^ (x e C[0, 1] I \mod{x, 6)\>e] ^0. (14) 

Then, for any reals < ti < t2 < ■ ■ ■ < tk < I, we show the convergence in 
distribution of the random vectors {X^^, X^, • • • , X^^) 

(x^, X^, • • • , X^) > (Xt^ , Xt2 , ■ ■ • , Xt^) . 

Dealing with finite-dimensional vectors, we can use the Cramer- Wold device 
if we designate and the characteristic functions of {Xl^,X^, ■ ■ ■ ,X^) 
and {Xti,Xt2, - ■ ■ i^t^) defined on W, we just have to show the point-wise 
convergence of the characteristic functions toward 

8.2. The Tightness of the Induced Family of Distributions 

We will confirm that the family of induced measure /i^ satisfy the tightness 
criterion l|14|) on C[0, 1]. First, define the random variables 

6„ = max ICn.fcl- 

0<fc<2"-i 

We apply the usual Borel-Cantelli lemma and state: there exists a set Q,' with 
P(ri') = 1 such that for every lo in fi', there is an integer n(a;) such that for all 
n > n(a;) we have 6n(ct;) < n. Let us set 

= {cj e r^' I n{Lu) < n} 

It clearly defines an increasing sequence of sets ^Iq C fl'i C ■ ■ ■ C 51'^ with 
limji^oo P(^^n) = 1- For any 77 > 0, there is a 71(77) such that for all n > n{r]), 
we have P(f^^) > 1 — 77/2. Then Considering the existence of the upper bound 
((8|). for every lu in with N > 71(77), can write 

n(i)) 

|Xf (C^)-Xf (C^)| < |*n,fc(i)-*n,fc(s)||en,fc(c^)| 

n=0 0<fe<2"(i)-i 
N 

+ Yl 271 •2-"*'.||.9||^. 11/11^ (15) 

n—n{rj)-\-l 
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The previous inequality actually holds for any iV > with the conventions of 
setting the elements ^n.k to zero when n > N. For all e > and all 77 > there 
exists an integer n(e, 77) > n{ri) such that 

00 

^ 2n.2-^.||5lL-||/IL<| 

n— n(e,?^) + l 

Now, for any given e > and 77 > 0, we chose a real M{e, rj) > large enough 
so that for any w in f2 

P( max 6„H>Af(e,77)) < ^ 

and we finally define the set fig^^ as 

^t,r, = < w e n) I max bn{uj) < M(e, 77) 

V ) /y n<n{e,ri) 

Every element Vl/n,fc is a continuous function defined on a compat support Sn.k 
in [0, 1]. As a result, the finite set of functions for 77, < n(e, 77) is uniformly 
equicontinuous: for any given e > and rj > 0, there is a real 6{e,ri) > such 
that 

Vt,se[0,l], \t-s\<d{e,r^) => 

V77, k with { I ^ ^ ^"ilV^) , l^^At) - ^nAs)\ < . 

For all e > and all 77 > 0, writing the inequality (fT5|) on fl^^jj yields to 
Vcj e (5 < (5(e, 77) ^ VAT e N, \mod{X^ (lu), 6) \ < e 

As we have defined i~le,ri so that P(0(:^^) > 1 — 77, we have shown that 

Ve,77>0, S<S{e,r]) sup P ( |77iod(A^(w), (5) I > e ) < 77 , 

which proves the tightness criterion on the Wiener space C[0,1] for the family 
of continuous processes A^. 



8.3. Point-wise Convergence of the Characteristic Functions 

For any (Ai,A2,-- - , Xr) in W, evaluating 'yf^, the characteristic function of 
(A,^,A,^,...,A,^), yields 
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We use the definition ([8|) of the partial sum in terms of the basis elements 
k to explicit the calculation of '■^^ on , the space of admissible functions 
for'x^: 



5 AT 



Remember that each coefficient ^n,k is independently distributed according to 
a normal law Af{0, 1) in the representation of . We have then 



'^(Ai,A2,---,A.)= n n / e^'-'=(S^'''*-''Mp(dC„,fc). 
Therefore, we can compute each terms of the previous product: 

JR V27r 



= exp - i(Ap^'„,fc(ip))^) . 
Back to the formulation of 'j^n , we end up with the analytical expression 

^^(Ai,A2,--- ,A,)=exp(--^^ XpX,p^{tp,s,)) , 

P 9 

where we have used the short notation for the covariance of X^ 

N 

n=0 0<fc<2"-i 

With the covariance calculation of section [71 we have demonstrated that for t, 
s in [0, 1] 

lim s) = E{Xt ■ Xs) = git)g{s) ■ h{min{t, s)) = p{t, s) . 

For any (Ai, A2, • • • , A^.) in R'^, it directly entails the point-wise convergence of 



r r 

lim 'g'^(Ai, Aa,- • ■ , Ar) exp ( - - V V ApA, p(tp, s 



7 AT 



P 9 
"^(Ai, A2, • • • , Ar) , 

where we remark that ^ is the characteristic function of the random vector 
{Xti,Xt2,- ■ ■ ,Xt^). 

By the Cramer- Wold device, the point-wise convergence of 'i^^ toward proves 
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the convergence in distribution of any finite-dimensional random vectors X^^, X^, ■ ■ ■ , 
toward the random vector Xt-^ , Xt2 , ■ ■ ■ , Xt^ 

{■^u.^ ^t^i ■ ■ ■ ' -^u) * {^ti , Xt2 , • • • 7 Xt^) . 

It also proves that X is an exact representation of the Gaussian Markov process 
X since, as the almost-sure path-wise limit of X^ when N tends to infinity, X 
follows the same law as the law of the process X. As we have already established 
the tightness of the family of induced distributions , it is enough to prove 
the weak convergence of /z^ toward /i on the Wiener space, which is equivalent 
to the convergence in distribution of the continuous process X^ towards X . 

Appendix A: Gaussian calculation 

We want to establish the analytical results given in section 13.31 Assuming < ty < ^s, we want 
to compute p{Xty — y \Xt^ — x^Xt^ — z). We use the Markov property as mentioned in the third 
section 

p(Xty = V \Xt^ = X, Xt, = z) = — — 

P{Xt^ = z \Xt^ = x) 

to get its analytical expression in terms of g and h 



(_y ( { ' 

\g(ty) s(tm)y 1 1 / Vs(tz) 



cxp 5 • — — , , oxp 



cxp 



with the expressions of the variance of Xt between any two consecutive points 



(16) 



<yl,y = h(ty)-h(t.) , <yl,. = Htz) - Hty) , <yl^,=h(t,)~h(t^) . (17) 
After factorization of the exponentials, the resulting exponent of expression is written 



y a; \ i / z y 



_L / 1 1 \ 2 1 / x_ 



+ 



(s(ti) 9(tx)) 



We then factorizc the term Cy so that we can write the exponent in the form 
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to obtain the canonical expression 



1 -yl 



git.) 



X + 



9(ty) 



— (— 



) ait.). 



2-g^(ty) 



g(*y) 
g{t.) 



git. 



2-g-'it.)ol^y 2-gHt,)al 



The quantity Q can be further simplify after being expanded 



2-g^it.) \'yl,,-'yl,, ' '^1-. 

1 <^x.y 2.1 1 

XZ - 



1 

XZ 



— ( — 

•^l,. \sit.) 



'^■'yl,.\ait.) git.), 

All terms except quadratic ones in y cancel out in the exponent, giving 



viXty =x,Xt^ =z) = 



I 



1 



/27r • gity-)- 



( S(ty} 
y 9(tx) 



■ x + 



(18) 



We sum up the expression il8l noticing it represents the distribution of a normal law J\f [x ^ity) , x'^ity 
whose parameters read 



XfJ-ity) = 



g(ty) 
git.) 



■ X + 



git. 



gjty) h(t.) ~ h{ty) 

git.) ' h(t.) - h(t^) 



X + 



gjty) hjty) - h{t.) 
git.) ' h(t.) - h(t^) 



xcrity) = g ity) ■ 



■ a ity) 



{h(ty) ~ h(t^))[h{ty) - h{t^)) 



h{t. 



Appendix B: Induced measures 

Finite-dimensional measures 

We first recall the definition of the 2^-dimcnsional vectorial space C" defined as 




We specify a given element x of by writting 
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a;(t) = ^ 5Z *i,fc(*)-a„,fc with a„,fceR, 

11 = 0<f!<2"-l 

and we remark that x can be viewed as a sample path of the process if we posit ^n.k(^) — f^Ti,k- 
We then introduce the following notations 



[n,k] = 2"'+k with <! [n'lj^a' 



r [0, 

1 [0, 



to enumerate indices between and 2 in the prefix order. With this convention, we introduce the 
positive reals — Iq < ti < ■ ■ ■ < t^j^ — 1 defined as 

t[n,k] ^ run^k ^ {2k + 1) 2~" for < n < and <k <n . 
For any collection of Borel sets Bq. B\. • • • , B^j^ in we define the cylinder sets 



which generate the cr-algebra B{C ). 

N 2^ 

We want to show that there exists a probability density p on IR such that we have for any sets 
{0} ^ Bo.,Bi.,--- ,B^N in 6(E) 



and that has the following simple expression in terms of the transition kernel p associated with 
the Gaussian Markov process X 

2^-1 

p^ {x-i,-- ■ , x^n) ^ p((a;A>fi, ifc+-i)|(a:fc, tfc)) . (19) 

k = Q 

We start noticing that, by construction of , we have 

P{CTi,fc e dar^^k) ^ 'P{Xm^^^ ~ p.{nirv,k) £ dX[n .k]\Xl ^ ^ X[r^-l,2k]. ^r^^^ ^ a: ^ 2 fcfl] ) 

with the usual definitions for ln,k, fn,k and rrtn^k and with dx^^j^-^ — (Tn,k ■ dan^k- As the Gaussian 
variables are all independent in the definition of X^ and remembering the definition of the 

probability density q in scction l3.3l we define for AT > as 

p^ {x-i;- ■ ■ , x^n) ^ p(^(x^N , ^^1,0) I (xq., ^i,o))<?((3:[o,i] , '^1,0) I (3^0, h,o), (X2N . '^1,0)) - ■ ■ 

2^-1-1 
fc = 

and we will show by recurrence on that p^ as the expected expression JlQI l. 

The basis statement is obvious for N — 0. As for the inductive step, let us write the probability 
density p^'^^ as 

N+l , . 

p ^ (a;i, ■ ■ ■ , a:;2iv+i ) ^ 

2«-l 

p^(a:i, ■ - - , aj^iv) 9 ((^^[iv+i,^] , "^iv+i,fc ) | Ca;[jv,2fc] : ^JV+i, A; ) ^ (a;[iv,2/^i] , '^iv+i,fc)) ■ 

By our recurrence hypothesis, we have then 
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JV+l/ N 



fc=0 fc=0 

But by definition of q, wc have 



p((2;[iV,2fcfl]. ^'JV+l,?;:) I (3;[Af,2fe] J ^iV+l,fe))) 

The product of the denominators in the previous expressions cancels out the term in so 
that we have proven that 

2^+1-1 

p^^'^{xi, - - - , 3:2^+1) ^ p((2^fcfi, ^fcfi) I ifc)) . 



Characteristic Functionals 

We want to express the characteristic functionals of the induced measures and of the Wiener 
measure fi. In view of this, we recall that, by the Riesz representation theorem, the dual space of 
C[0, 1] is the space M[0, 1] of finite measures on [0, 1]. For x in C[0, 1] and z-* in M[0. 1], we write 
the duality product 

{x, ly) ^ I x{t) du{t) . 





Considering first the induced measure /^jv, the characteristic functional 'i^jv is defined on -/\/[0, 1] 
the Fourier transform of p.j^ on C[0, 1] 

^jv(!^) - / e'^^' d^^(x) with ueM[{),l]. 
J c 

By construction of the process and independence of the random variables ^n^h-, we have 
The random variables ^ri,k being Gaussian of law A/^(0, 1), we furthermore have 

This allows to finally write the expression of the characteristic functional S^jv 



Ti = 0</s<2^i-l 

If we express the product of duality, we can formulate the functional "^iv in its common form 
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exp 

'0 J Q 

where we made apparent , the correlation function of the process 



Ti = 0<A]<2"-1 

From the convergence results of our expansion, we directly have that the characteristic functional 
of the Wiener measure f 



(20) 



"^(i/) = cxp ^ - i y J p(t,s) du{t)du{s)J , 

with p the continuous correlation function of the Gaussian Marlcov process X : 

p(t,s) = E{Xt ■ Xs) = g{t)g(s) ■ /i ( min(t, s)) . 
Now let ti < ^2 < • • ■ < tfc be some reals in [0, 1], we posit the measure u as 

k 

^ = ^ K^tp , 
p=i 

where St^ denotes the Dirac distribution concentrated at tp. If we inject the expression of i/ in the 
result t20t . we find the expression of the characteristic function of Xt-^ , Xt2 , ' ' ' , ^tf^ as one might 
expect. 

Then, assume the distribution f admits a density 9 in L^(0, 1) with respect to the Lebesgues 
measure, we have 



J 

Thanks to the auxiliary orthonormal basis ^n^k: we can further write 



dt , 



,k 



which directly leads to the following simple expression for the characteristic functional 



i/)=exp(^-i f (/(t) f gis)e{s)ds) dt^ 



References 

[1] BiLLiNGSLEY, P. (1999). Convergence of probability measures, Second 
ed. Wiley Series in Probability and Statistics: Probability and Statistics. 
John Wiley & Sons Inc., New York. A Wiley-Interscience Publication. 
pRMR1700749 (2000e:60008) 



Thibaud Taillefumier/ A Discrete Construction for Gaussian Markov Processes 26 



[2] Cameron, R. H. and Donsker, M. D. (1959). Inversion formulae for 
characteristic functionals of stochastic processes. Ann. of Math. (2) 69, 15- 
36. MRMR0108850 (21 #7562) 

[3] CiESiELSKi, Z. (1961). Holder conditions for realizations of Gaussian pro- 
cesses. Trans. Amer. Math. Soc. 99, 403-413. MRMR013259]_(24#A2431)] 

[4] DoOB, J. L. (1949). Heuristic approach to the Kolmogorov-Smirnov theo- 
rems. Ann. Math. Statistics 20, 393-403. MRMR0030732 (11,43a) 

[5] HiDA, T. (1960/1961). Canonical representations of Gaussian processes and 
their apphcations. Mem. Coll. Sci. Univ. Kyoto. Ser. A. Math. 33, 109-155. 
MRMR0119246 (22 #10012) 

[6] HOTELLING, H. (1933). Analysis of a Complex of Statistical Variables into 
Principal Components. Journal of Educational Psychology 24, 7, 498-520. 

[7] Kac, M. AND SiEGERT, A. J. F. (1947). An explicit representation 
of a stationary Gaussian process. Ann. Math. Statistics 18, 438-442. 
'MRMR0021672 (9,97a) 

[8] Karatzas, I. AND Shreve, S. E. (1991). Brownian motion and stochastic 
calculus, Second ed. Graduate Texts in Mathematics, Vol. 113. Springer- 
Verlag, New York. MRMR1121940 (92h:60127) 

[9] Karhunen, K. (1946). Zur Spektraltheorie stochastischer Prozesse. 
Ann. Acad. Sci. Fennicae. Ser. A. I. Math.-Phys. 1946, 34, 7. 
MRMR0023012 (9,292h) 

[10] Levy, P. (1948). Processus Stochastiques et Mouvement Brownien. Suivi 
d'une note de M. Loeve. Gauthier-Villars, Paris. MRMR0029120 (10,551a)J 

[11] Loeve, M. (1963). Probability theory. Third edition. D. Van Nostrand 
Co., Inc., Princeton, N.J.-Toronto, Ont.-London. MRMR0203748 (34 #3596)] 

[12] LUMLEY, J. L. (1967). The Structure of Inhomogeneous Turbulent Flows. 
Atmospheric turbulence and radio propagation, 166-178. 

[13] Rogers, L. C. G. and Williams, D. (2000). Diffusions, Markov pro- 
cesses, and martingales. Vol. 1. Cambridge Mathematical Library. Cambridge 
University Press, Cambridge. Foundations, Reprint of the second (1994) edi- 
tion. MRMR1796539 (2001g:60188) 

[14] Taillefumier T., M. M. A haar-like construction for the ornstein uhlen- 
beck process. Accepted in J. Stat. Phys.. 



